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The fragmentation of expanding shells II: 
Thickness matters. 
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ABSTRACT 

We study analytically the development of gravitational instability in an expanding 
shell having finite thickness. We consider three models for the radial density profile of 
the shell: (i) an analytic uniform-density model, (ii) a semi-analytic model obtained 
by numerical solution of the hydrostatic equilibrium equation, and (iii) a 3D hydro- 
dynamic simulation. We show that all three profiles are in close agreement, and this 
allows us to use the first model to describe fragments in the radial direction of the 
shell. We then use non-linear equations describing the time-evolution of a uniform 
oblate spheroid to derive the growth rates of shell fragments having different sizes. 
This yields a dispersion relation which depends on the shell thickness, and hence on 
the pressure confining the shell. We compare this dispersion relation with the disper- 
sion relation obtained using the standard thin-shell analysis, and show that, if the 
confining pressure is low, only large fragments are unstable. On the other hand, if the 
confining pressure is high, fragments smaller than predicted by the thin-shell analysis 
become unstable. Finally, we compare the new dispersion relation with the results of 
3D hydrodynamic simulations, and show that the two are in good agreement. 
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1 INTRODUCTION 

Shells and bubbles are common morphological features 
in the interstellar media of galaxies, and have been 
many diff er ent w aveleng ths. 



observed at many diff er ent w aveleng ths. For instance, 
iMcClure-Griffiths et all l|2002l ') and lEhlerova fc Palou^ 
( 2005|) have identified ~ 1000 expanding shells in HI surveys 
of the Milky Way. HI shells have also been found in the LMC 
iJKimet al.lll998l ). the SMC ( Hatzidimitriou et all l2005ll. 
and other ne arby galaxies. In the infrared, IChurchwell et al.l 
l|2006l . bOOTl ) have assembled a catalogue of ~ 600 shells 
found in the GLIMPSE survey of our Galaxy, using the 
SPITZER space telescope. Shells have been detected in th e 
Wisconsin Ha mapper data fWHAM: iHaffner et aLll2003h . 
Shells have also been studied in millimetre molecular line 
emission, radio continuum and X-rays. 

It is generally believed that a substantial fraction of 
these shells is formed by feedbac k from massive stars . In ad - 
dition, it has been suggested by lElmeereen fc Ladal (jlOTTI) 
that expanding shells can trigger star formation by the 
collect- and- collapse mechanism. In this mechanism, a mas- 
sive star (or a group of stars) injects energy into the inter- 
stellar medium in the form of stellar winds and ionising ra- 
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diation, and creates a bubble of gas with high temperature. 
The bubble expands due to its high internal pressure, and 
sweeps up the ambient medium into a dense shell. The shell 
cools down, becomes gravitationally unstable, fragments and 
forms a new generation of stars. If the new generation in- 
cludes massive stars, the whole process may repeat itself, 
leading to sequential self-propagating star formation. Vari- 
ous modific ations to the collec t-and-collapse mechanism are 
reviewed bv lElmegreerJ (119981 ). 

The collect-and-collapse mecha nism has been tested ob- 
serva t ionally in a s eries o f pap ers bvlDeharveng et al.l (120031 . 
l2005l . l2006l . l2008l . |2009| ) and IZavagno et al.l l|2006l ). Using 
data at various wavelengths, they identify and study sev- 
eral objects in which star formation a ppears to hav e been 
triggered at the borders of HII regions. IWatson et al.l (|2008l . 
2009) analyse six shells from the GLIMPSE survey, and de- 
termine the properties of young stellar objects concentrated 
in these shells. 

The gravitational instability of an expanding infinites- 
imally t hin shell h as be en studied theoretically by several 
authors. IVishniad l|l98a ) derives a dispersion relation (per- 
turbation growth rate as a function of wavenumber) for an 
infinitesimally thin shell, using decomposition into spherical 
harmonics. A ver y similar dispersion relation is obtained by 
lElmegreenl l|l994 ) using perturbation analysis of linearised 
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two-d imensional hydrodynamic equations. IWhitworth et alj 
([199J) use the equation of motion for a two-dimensional 
fragment forming on the surface of a shell, and find expres- 
sions for the size of the first and most unstable fragment, and 
for the time at which the instability starts. All these analy- 
ses use a very similar physical model (an infinitesimally thin 
shell which expands into a uniform medium), and although 
the mathematical description in each of them is very differ- 
ent, the results regarding fragment sizes and growth rates 
are very similar. 

We have tested these analytic predictions for the growth 
rate of gravitat ional instability, based on the thin-shell ap- 
proximation, in lDale et al.l (|2009| . hereafter Paper I). We use 
two very different three-dimensional hydrodynamic codes 
(an Eulerian AMR code and a Lagrangian SPH code) to sim- 
ulate the evolution and subsequent fragmentation of expand- 
ing self-gravitating shells. Our setup differs slightly from the 
one studied by Vishniac, Elmegreen and Whitworth, in that 
our shell does not accrete mass from the ambient medium. 
This modification requires only minor changes to the linear 
theory of the thin shell gravitational instability, which we 
describe in i)3.1l 

We make this modification because an accreting shell i s 
prone to the Vishniac dynamical instability (| VishniadI 1983h . 
which would grow quickly and make an analysis of gravita- 
tional instability impossible. On the other hand, we want to 
have the option to keep the shell thin, in order to test the 
thin-shell approximation. Therefore, we fill the shell inte- 
rior and exterior with a hot rarefied gas and the pressure of 
this gas confines the shell. However, because of its low den- 
sity, the ram pressure of the accreted gas, which is crucial 
for the development of the Vishniac instability, is negligi- 
ble. We also keep the pressure of the rarefied gas constant 
throughout the evolution, so that the shell is effectively in 
free fall, and development of the Ray leigh- Taylor instability 
is suppressed. 

In Paper I, we find excellent agreement between the two 
numerical codes. However, for simulations with low confin- 
ing pressure, in which the shell becomes thick, the simula- 
tions do not agree well with the predictions of the thin-shell 
approximation. The agreement is better for simulations in 
which the confining pressure is chosen so that the shell thick- 
ness stays approximately constant during its evolution. 

In this paper (Paper II) we study the gravitational in- 
stability of a thick shell analytically, in order to understand 
how a thick shell fragments, and why the fragmentation dif- 
fers from the predictions of the thin-shell approximation. 
We show that the confining pressure, which defines the shell 
thickness, is an important factor in determining the range 
of unstable wavelengths. Therefore, we call the new model 
pressure assisted gravitational instability (PAGI). 

The outline of the paper is as follows. In ^we derive a 
semi-analytic description of the equilibrium radial profile of 
the shell, and compare it with our 3D numerical simulations, 
and with the simple uniform-density model of the shell that 
we adopt in the subsequent analysis. Section ij3] deals in 
detail with modelling gravitational instability in a shell. We 
briefiy derive the dispersion relation for a non-acreting shell, 
using the thin-shell approximation; we give a description 
of a fragment forming in a shell with non-zero thickness, 
and compare it with the description used in the thin-shell 
approximation; finally we derive a new dispersion relation 



for a thick pressure-confined shell. In ^ we compare the 
new dispersion relation with growth rates measured in our 
3D numerical simulations, in Sj5]we discuss the astrophysical 
consequences of our results, and in 
conclusions. 



we summarise our 



2 RADIAL DENSITY PROFILE OF THE 
ISOTHERMAL SHELL 

Since we wish to study the infiuence of the shell's thickness 
on its gravitational instability, it is essential first to under- 
stand the shell's equilibrium radial structure. In Paper I we 
set the initial radial density profile of the shell using a simple 
local model for the hydrostatic equilibrium of the gas layer, 
and assuming that the shell radius is large compared with 
its thickness: 



p{R) — Po sech 



2ttGpc 
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(1) 



Here R is the radial coordinate, p^ is the maximum density, 
which appears at the radius R^ , and Cg is the isothermal 
sound-speed. 

A more realistic density profile can be obtained by drop- 
ping the assumption that the density profile is symmetric 
about _Rq and solving the equation of hydrostatic equilib- 



rium in the form 
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Here M{R) is the mass within radius R, and q is a dimen- 
sionless number between and 1 which has to be found 
numerically (see below). The first term on the RHS of Eqn. 
[2]represents the net gravitational acceleration at radius i^. 
The second term on the RHS represents the gravitational ac- 
celeration of the material at radius i^ , which is the radius 
that contains a fraction a of the total mass of the shell. This 
gravitational acceleration is subtracted from the net gravi- 
tational acceleration, on the grounds that it contributes to 
the bulk deceleration of the centre of mass of the expanding 
shell, rather than its radial profile. 

Switching to a Lagrangian formalism (i.e. using Af as 
the independent variable) Eqn. [2] becomes 
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where ^' = d^/dp and {" = d^(,/dp'^. 

The shell density profile can be found by solving Eqn. 
((6]) numerically, but we also need to determine the appro- 
priate boundary conditions and the value of ct. The shell 
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is confined by the external pressure P^-^^ on both its in- 
ner and outer surfaces, and therefore the density at these 
surfaces has to be p^ — P^xt /^s • Inserting this into the 
mass conservation law we obtain the inner boundary con- 
dition in the form (C'C^) = A^TOTCs/(4vr-Rf Pext')- The 
value ^(/^ = 0) is unknown, but it can be found by iterating 
until the boundary condition is also fulfilled at the outer sur- 



face, 



{ee 



= M^a^cl/{'i%RlP^ 



The parameter 



a can be found by adding another level of iteration in which 
a is varied until it fulfils the initial assumption ^(/i = a) = 1. 
Figure [T] compares the radial density profile of the 
shell in the AMR run from Paper I with P^xt = 
10~^^ dynecm"^, at time 7.2 Myr, with the semi-analytic 
profile obtained using Eqn. ((6]) and the same shell param- 
eters. Both profiles agree very well with each other, which 
means that the assumption of hydrostatic equilibrium is cor- 
rect and that the shell is well resolved in the numerical simu- 
lation. The time evolution of the shell radius and its FWHM 
are shown in Figure [21 and this also demonstrates excellent 
agreement between the numerical and semi-analytic treat- 
ments. In ^ we assume that the shell has uniform density 
with half-thickness given by 



Zj^C^ 



2-Pf,xt + ■^G'E^ 



(7) 



where Eq = M.j.q^/{A-kR^) is the shell surface density. 
This equation is obtained by splitting the shell into two 
identical layers and then equating the external pressure 
force and gravitational force acting on one layer due to 
the other one with the internal pressure force in between 
the layers. Figures [l] and [2] also display the density pro- 
file, mean radius and FWHM of this uniform-density shell 
model, and show that it is a good approximation to the 
more detailed models. The agreement between the numeri- 
cal, semi-analytic and uniform-density models is even better 
for runs with higher external pressures P^xt = 10~^^ and 
P^^^ = 5 X 10"" dynecm"^. 

Figure [2] also shows the evolution of the parameter a. 
It is close to 0.5 in the beginning of the expansion when the 
shell is thin and its profile is symmetric, and only decreases 
by ~ 10% subseque ntly. This result suppor t s the assumption 
of a = 0.5 made bv lWhitworth fc Francis l|2002l '). 



3 GROWTH RATE OF THE GRAVITATIONAL 
INSTABILITY 

The main aim of this paper is to derive the dispersion rela- 
tion for gravitational instability of a thick shell, i.e. pertur- 
bation growth rate as a function of perturbation wavenum- 
ber. For this purpose we model a forming fragment as a 
uniform-density spheroid whose radial excursions are de- 
scribed by non-linear differential equations. i|3.1l brieflv sum- 
marises the derivation of the thin-shell dispersion relation 
for the case of a non-accreting shell. t|3.2l describes the 
spheroidal model for a fragment forming in a thick shell, 
and explains in detail the differences between this model 
and the model underlying the thin-shell approximation. tj3.3l 
presents the derivation of the PAGI dispersion relation and 
compares it with the thin-shell dispersion relation. 
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Figure 1. Radial profile of the shell at t = 7.2 Myr in 3D AMR 
simulation (plus symbols), semi-analytical model (solid line) and 
uniform shell with the same column density (dashed line). 




Figure 2. Evolution of the shell FWHM. The dash-dotted line 
denotes the position of the shell peak density, the other lines 
below and above it show positions of the the shell inner and outer 
edge measured at half of the peak density for the three models: 
3D AMR simulation (plus symbols), semi-analytical model (thin 
solid line) and uniform shell (thick dashed line). The thin dashed 
line shows the evolution of the parameter a. 



3.1 Thin shell approximation 

Consider an infinitesimally thin shell with radius R and sur- 
face density E, which expands with velocity V into a vac- 
uum. The shell is isothermal, with sound speed Cg . The equa- 
tions describing transverse flows of gas inside the shell arqj. 



9E 
dt 



+ PEVt 
,30. 



V 
+ 2^- =0, 



E7?V = -cf VtE - EV$ - EVn , 



dt 



4TvGT.S{r ~ R) . 



(8) 

(9) 
(10) 



Here O is the angular velocity of transverse flows, $ is the 
gravitational potential, and V^, and V^,- denote the two- 



1 In lElmegreenI |T99i) and IWiJnsch fc Paloud lHooJ) the shell 
was assumed to accrete mass from the ambient uniform medium, 
and this resulted in a factor of 4 in front of the last term in the 
equation of motion, Eqn. |9] 
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dimensional gradient and divergence, respectively. The in- 
stability growth rate is obtained by substituting plane- wave 
perturbations parametrised by the dimensionle ss wavcnum- 
ber I = kR (see [E lmcgrcei] Il994 IWiinsch fc Paloua 200 li 
and Paper I for details), 



E = Eq + E; cos(;6l) , 

n ^ fli sm{W) , 

$ = <&o + *i cos(Z6') . 



(11) 



Here Eq and "I>q are unperturbed values of the surface den- 
sity and gravitational potential, Ej, Qi and "I>; are the per- 
turbation amplitudes, and 9 is an angular coordinate on 
the surface of the shell in a direction perpendicular to the 
plane wave. Substituting these perturbations (equation lllf) 
into the linearised equations (|SHTDJ, and integrating equa- 
tion pup over the thickness of the shell, gives a set of linear 
differential equations for the amplitudes E; and SI; 



di 



E, 



2V 
R ' 



2ttG 
R ' 
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R 




(12) 



The set of equations (I12|l has two eigenvalues 



,(1,2) 
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(0 = -^±\ Tm + 



V^ 2ttGT.„1 c. 



2R \l AR'^ R 

and two related eigenvectors 



o(l>2) 



,(1,2) , 2V 



"w (1^) 



(14) 



Any harmonic plane-wave perturbation of the shell de- 
scribed by amplitudes E; and Qi can be written in the form 



Qi 



= ^er+Sef), 



(15) 



where A and B are real numbers. Its evolution in the linear 
regime is then given by 



^ , = Aep' exp(a;(''i) -|- Bef^ exp(a;(''i) . 



(16) 



If a perturbation is given by the eigenvector (2) only (i.e. A 
— 0), the mode is always damped (since cj'^' is always neg- 
ative). On the other hand, if some perturbation is described 
by eigen-vector (1), it will grow if the gravity term in oj'^' 
is large enough. 

3.2 Fragment in the thick shell 

It turns out that the thin-shell approximation is a good de- 
scription of the linear instability at wavelengths larger than 
the shell thickness, provided the shell thickness is uniform. 
However, in Paper I we find that the shell thickness varies 
with the surface density of the perturbation, even if its am- 
plitude is relatively small. The spatially varying shell thick- 
ness has two effects, (i) In a low-pressure environment, the 
gas at the centre of a fragment expands in the direction per- 
pendicular to the shell surface, and this results in a slower 
growth rate and a smaller range of unstable wavelengths, 
compared with the predictions of the thin-shell approxima- 
tion, (ii) In a high-pressure environment, the external pres- 
sure helps to compress fragments from the side, and this re- 
sults in a higher growth rate and a larger range of unstable 



wavelengths, compared with the predictions of the thin-shell 
approximation. 

In order to take the variable thickness of the shell into 
account, we model a fragment as a uniform-density oblate 
spheroid with major axis r and minor axis z, embedded in an 
ambient medium with pressure ^ext- Radial excursions of 
the spheroid are then controlled by the equations of motion 



r{r, z) 



and 



'i{r,z) 



3Gm j r cos ^{z/r) 



z/r 



(r^ — 2^) 



2^3/2 



(r2 



20tv Pj^^^rz 
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5cf 



(17) 



-3Gm 



(z/r) 



(r2 



(r2 
207rP„- 



^2)3/2 

rf 5c; 



-^ (18) 



3m z 

(see iBovd fc Whitworthll2005l ) . 

We assume that the mass of the spheroid, m, is con- 
stant, and originates from a comoving circular patch on the 
shell, with constant angular wavenumber I, and hence radius 
nR(t)/l. 

The evolution of the spheroid can be followed 
by solving equations (|17l) and (|18|) numerically, as in 
I Boyd fc WhitworthI (|2005|). However, in order to derive an 
analytic dispersion relation for PAGI, we define i = as 
the time when fragmentation starts, and then we assume 
that for a certain time, te, thereafter the spheroid collapses 
with constant acceleration r^ given by equation (|17|l with 
r = Tq. Here r^ — ttR/I is the initial radius of the fragment, 
and its initial expansion speed is Vg — nV/l, where V is 
the radial velocity of the shell. The subscript e denotes the 
fraction by which the fragment radius shrinks during t^ (i.e. 
r(ie) = (1 — e)rg). The evolution of the spheroid radius is 
then given by 



r{t) ^rg + rgt+-rgt^ . 
and hence t^ is given by 
rii? + 2rgU + 2erg = . 



(19) 



(20) 



In the thin-shell analysis, the instability growth rate is 
defined as the inverse of the time it takes for the perturba- 
tion described by a given eigenvector to grow by a factor e. 
Since eigenvectors are not available for PAGI, we define the 
instability growth rate simply as 

a;, = i . (21) 

Since u^ depends on the choice of e, its absolute value is 
somewhat arbitrary. However, relative values of lu^ for dif- 
ferent wavenumbers are well defined. Moreover, the range of 
unstable wavenumbers is independent of e and can therefore 
be compared directly with the range predicted by the thin- 
shell dispersion relation a;rpjjj;.j. Combining Egns. l^ and 1201 
we obtain the PAGI dispersion relation in the form 



(22) 



Figure [3] illustrates the models used in the two analy- 
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by eigenvectors (Egn. [T4)) . In the PAGI analysis, the frag- 
ment starts life as a uniform-density spheroid expanding 
with the shell, (iii) The self-gravity term reflects the vertical 
height of the spheroid, 2^. For a spherical fragment {z^ — Vg), 
the term in the square bracket is 2/3r„. As the fragment 
becomes flatter, this term increases, and in the limit of an 
infinitesimally thin fragment {zg = 0) it becomes Tr/2rg (an 
increase of 3tt/4 ~ 2.36). (iv) The fourth term in Egn. I23l 
(line three) represents the effect of external pressure, and 
therefore has no equivalent in Eqn. 1131 



Figure 3. Comparison of fragment models in the thin shell ap- 
proach (right) and the thick shell approach (left). The approxima- 
tions are: sine-shaped surface density and velocity perturbation 
(thin) vs. uniform oblate spheroid (thick); and solution of lin- 
earised hydrodynamic equations (thin) vs. constant acceleration 
coming from non-linear equations of motion (thick). 



ses of shell gravitational instability. In the thin-shell anal- 
ysis, the fragment is modelled as a sine-shaped perturba- 
tion of the surface-density and the velocity, on an infinites- 
imally thin shell. Its growth rate is obtained by solving lin- 
earised hydrodynamic equations. In the PAGI analysis, the 
fragment is modelled as a uniform-density oblate spheroid 
contained by a constant ambient pressure. It is assumed to 
contract with constant acceleration given by the non-linear 
equations of motion evaluated at the onset of instability. 



Due to the dependence of the absolute value of the 
PAGI dispersion relation on e it is not possible to compare 
ujc to (j-'thin directly (we remind the reader that only the 
absolute value of Ue depends on e, the range of unstable 
wavelengths as well as the relative growth rates of individual 
modes are e-independent). However, we can deflne a critical 
pressure P-j-nm such that with P^xt ~ ^thin the range of 
unstable wavenumbers predicted by the PAGI analysis is the 
same as that predicted by the thin-shell analysis. Then, for 
this choice of -Pexti ^^ ^^'^ ^ ^^^ which Eqns. (I23|l and p3|l 
predict very similar growth rates {u}f{l,P^-^jj^) ~ uj^-^jj^{l)). 
This is shown by Figure [J] for the shell from Paper I in two 
states: the initial radius _R = 10 pc and the maximum ra- 
dius -R = 23 pc. It can be seen that ^^{1, Pthin) — '^thin(0 
if we set e — 0.1, for both shell radii. Moreover, this value 
of e results in both dispersion relations being close to each 
other throughout the whole shell expansion and we therefore 
adopt this as our default value of e. 



3.3 Dispersion relation of the thick shell 

The final form of the dispersion relation is obtained by sub- 
stituting rr\, = ■''"'"o^^i '"0 — 1"-^/^! ^0 = '^K/^ (where I^, i^ 
and Vg are - respectively - the surface-density, radius, and 
radial velocity of the shell when the instability commences) , 
and r„ = r{r„, zQ (Eqn. [13 with 2^ = Ec^ / {2P^^^ + ttGE^) 
(Eqn. [3. This yields 



2eHg \\2eH„J 
37rGH, 



4e 



Vg cos ^ (^/r: 



% 



+ 



- (''0 - ^^) 



^^Pext ^s 



2^3/2 



37r2ei^2(2P^^^ 



-^GE2) 



2Tv'2eB^ 



(23) 



Equation (|23|) has a very similar structure to the thin- 
shell dispersion relation (Eqn. [13]), in that the first two terms 
on the right hand side (line one) represent stretching due to 
expansion, the third term (line two) represents self-gravity, 
and the final term (line four) represents internal pressure. 
However, there are significant differences, (i) The terms all 
depend on the value of e. We revisit this issue below, (ii) 
The stretching terms have different numerical coefficients 
because the initial conditions for the analyses are different. 
In the thin-shell analysis, the initial perturbation consists of 
harmonic waves in both surface-density and velocity, given 



Figure [5] compares the PAGI dispersion relation for 
different values of Pj-xt with the thin-shell dispersion rela- 
tion. We have used parameters for a shell having a mass of 
2 X 10'* Mq, (a) with radius 10 pc and outward radial veloc- 
ity 2.2kms~*, and (b) with radius 23 pc and zero outward 
radial velocity; this is the same as the shell that we simu- 
lated in Paper I. Increasing P^xt results in a larger range of 
unstable wavenumbers, and faster growth rates, but even as 

MAX 

^EXT ~^ °°i the maximum unstable wavenumber, (p^(,j re- 

MAX MAX 

mains finite. Specifically, for Pj,xt = 0' 'pagi ~ ^-''^''^thin ' 

MAX MAX MAX 

and for Pext = 00, /p^^^j = 2.21^^^^^^ , where l^^^,^ is the max- 
imum unstable wavenumber according to the thin-shell dis- 
persion relation; evidently an infinitesimally thin shell (i.e. 
= 00) does not give the same dispersion 



PAGI with Pex 

relation as the thin-shell model. 



__ MAX MAX 

Figure [d] shows how the ratio 'pagi/^thin varies with 
PgxT- The vertical line denotes the critical external pres- 
sure, PcHiT = ''^GT^/2, at which the derivative of the shell 
thickness (z^g) with respect to its surface density (Eg) is 
zero, i.e. the configuration in which external pressure and 
self gravity make equal contributions to containment of the 
shell, and the shell thickness is only weakly dependent on 
its surface-density. This is why Pcpu^ is approximately equal 

to Ppjjji,, (where ^pagi/^thin = !)■ With this external pres- 
sure, as a fragment starts to form, and its surface-density 
increases, its thickness does not change much, and therefore 
the external pressure force remains largely perpendicular to 
the surface of the shell and does not make a significant con- 
tribution to the lateral squashing of the fragment. 
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Figure 5. Perturbation growtli rate for diflferent values of external pressure for two states of the shell from paper I: R 



V = 2.2 kms ^ (left), i? = 23 pc, V = (right). Alternating dashed and solid lines denote u),, for P^ 



10 



-10 



: 0, 10" 



10" 



10" 



10 pc, 
10-11, 



and 1 dyne cm ^ (in this order for growing range of unstable wavenumbers) . The dash-dotted line shows perturbation growth rate 



of the thin shell a;„ 



The horizontal dashed line at u; = represents the separation between unstable and damping modes. 



B. = 10 pc 

F„„., = 1.37 X 10-'^ dyne cm-^- 




Figure 4. Comparison of the thin shell dispersion relation 
(Eq. I I13I I. solid lines) and the PAGI dispersion relation (Eg. I23I I 
for three values of e: e = 0.05 (dashed), e = 0.1 (dash-dotted) 
and e = 0.2 (dotted). Two states of the shell with parameters 
from Paper I are shown: the initial state {R = 10 pc) and the 
state with the maximum radius ij = 23 pc. For each state, the 
external pressure P^xt ~ ^thin ™^® calculated so that ui^-^-^-^il) 
and u)e(l, ^thin) S'^*^ *'^'^ same range of unstable modes. 



4 COMPARISON WITH NUMERICAL 
SIMULATIONS 

In Paper I we have studied numerically the evolution of a 
shell with mass A'/ = 2 x IO^Mq, temperature T = 10 K, 
initial radius i?ini = 10 pc and initial expansion velocity 
Vn-ii — 2.2 kms-1, embedded in a rarefied medium with pres- 
sure -Pext- After about ISMyr, the expansion stalls, at ra- 
dius 7?max = 23 pc, and the shell starts to collapse. 

We have run simulations with three different exter- 



nal pressures -P^xt ~ 10 



= 10" 



and P„ 



5 X 10^ dyne cm" ; the first one is as close to zero as the 
hydrodynamic solver permits, the second one is chosen so 
that the shell thickness varies as little as possible during the 
shell evolution, and the third one is the highest pressure for 




Figure 6. Ratio of the shortest unstable modes given by the thick 
shell and thin shell models as a function of the external pressure 
^EXT ■ '^^'^ vertical dotted line denotes value of Pq^i,^ for which 
the shell has the maximum thickness for a given surface density. 



which the shell thickness can be resolved with the available 
computational power. 

Figures [7] and |8] summarize results from the hydrody- 
namic simulations. In the low-pressure run, we find that the 
perturbation growth rates diverge significantly from the pre- 
dictions of the thin-shell dispersion relation cjrpjjj;,j(Z) - only 
the low wavenumbers grow. The agreement is better for the 
medium-pressure simulation. The high pressure run again 
diverges from the predictions given by tj^^j;,) [1) ~ the range 
of unstable modes includes higher wavenumbers than pre- 
dicted. Differences in relative growth rates among the three 
simulations result in different fragment sizes as can be seen 
in figure [S] 

The dispersion relation obtained from numerical sim- 
ulations exhibits several peaks in some cases (see the top 
middle panel of Figure [7|. Among them, only the first peak 
(at lowest wavenumbers) is suitable for analysis of mode 
growth rates. The other peaks represent higher harmonics 
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Figure 9. Surface density profile of the fragment from the sim- 
ulation with i^EXT ~ 10"^'^ dynecm"^ at i = 15 Myr (corre- 
sponding to the top middle panel of Figure [7] and the left panel 
of Figure [Sj. The profile is taken along the Hammer projection 
coordinate at a; = —0.7. The thick solid line is the actual sur- 
face density in the simulation, the other lines are a sequancc of 
its spherical harmonics representation with modes up to a certain 
wavenumber: I < 15 (dotted), I < 30 (thin solid), Z < 45 (dashed), 
I < 60 (dash-dotted) and Z < 95 (dash-doubledotted). 



of the first peak and their growth is induced by the growth 
of modes from the first peak. The higher harmonic peaks 
describe fragment shapes which are typically more concen- 
trated than pure spherical harmonics with wavelengths cor- 
responding to fragment sizes. This is illustrated by Figure [9] 
which shows the surface density profile across an arbitrary 
fragment in the simulation with P^xt = 10~^^ dynecm"'^ 
and a sequence of its representation by spherical harmon- 
ics which take into account only modes up to a certain 
wavenumber. It can be seen that the modes from the first 
peak only (Z < 30) are enough to obtain approximately the 
correct fragment wavelength. If the second peak (Z < 45) 
is taken into account the representation of the fragment is 
already very good. Taking into account modes with Z > 45 
results in further small improvements in representing the 
detailed shape of the fragment. 

The perturbation growth rates from the numerical sim- 
ulations agree much more closely with the predictions of the 
PAGI dispersion relation (Egn. [23|l . The agreement is best 
at intermediate times. At earlier times the numerical sim- 
ulation has not yet had sufficient time to relax and so the 
modes are poorly defined. At later times the gravitational in- 
stability has become strongly non-linear. However, the PAGI 
dispersion relation does tend to give a slightly larger range 
of unstable wavenumbers than the numerical simulations. 



5 DISCUSSION 

The PAGI dispersion relation predicts that the fragment 
mass function depends on the external pressure confining 
the shell. Shells expanding into a low pressure environment, 
for instance shells at higher galactic latitudes, should frag- 
ment preferentially on long wavelengths. This may result 
in a top-heavy core mass function (CMF). However, it re- 
mains an open question whether the massive fragments will 



form massive stars or whether they will continue to fragment 
into smaller pieces. In the low pressure simulation presented 
in this work, the large fragments include relatively dense 
cores (whose growth is represented by distinctive higher 
harmonics in measured mode growth rates), and the evo- 
lution exhibits rather merging of fragment cores than frag- 
ment splitting. However, this can be due to our physical 
model being too simple (isothermal equation of state, rel- 
atively small perturbations of initial density and velocity, 
no radiation, stellar feedback, etc.). To answer this question 
properly, more realistic simulations would be necessary. 

If a top-heavy CMF results in a generation of stars 
with a top-heavy IMF, this should be favourable to the 
production of further shell s, and hence to sequenti al self- 
propagating star formation (jElmegreen fc Ladalll977l ). How- 
ever, it is a question if such secondary shells expanding 
into low density gas are able to collect enough mass, cool 
down and fragment. An example of such a shell in a low 
density environment might be the Carina Flare super-shell 
GSH 287-1-04-1 7, which extends up to ~ 45 pc above the 
galactic plane (jFukui et a l. 1999; Dawson et al. 2008). 

In Paper I and in this paper, we have studied shells ex- 
panding into a hot rarefied gas. Such shells are confined only 
by their self gravity and by external thermal pressure. This 
model is applicable to shells which break out of a molecular 
cloud and expand into a rarefied warm intercloud medium. 
Our reason for focusing on this model is that we want to 
study pure gravitational instability. By reducing the den- 
sity of the ambient gas and hence making its ram pressure 
negligible, we suppress the typically much faster Vishniac 
instability. 

For a shell expanding into the relatively dense gas of 
a molecular cloud, ram pressure plays an important role in 
compressing the outer surface of the shell. The Vishniac in- 
stability is then inevitable, because ram pressure acts purely 
parallel to the direction of shell expansion, whereas ther- 
mal pressure acts perpendicular to the surface of the shell. 
Furthermore, ram pressure cannot compress a fragment lat- 
erally, and we therefore expect that in this situation the 
growth of small fragments would be slower than predicted 
by the PAGI dispersion relation (Eqn. [23]). 

The PAGI dispersion relation predicts a range of unsta- 
ble wavenumbers which is systematically larger by ~ 25% 
than the range obtained from numerical simulations. This 
discrepancy appears at high wavenumbers, so it is unlikely to 
be due to the curvature of the shell (which is not taken into 
account in the derivation oi lu^), because the high wavenum- 
bers are least affected by it. The most probable source of this 
discrepancy is the approximation of spheroid evolution by 
collapse (or expansion) with constant a cceleration (Eq. I19II. 
This is certainly not true, as shown bv lBovd fc WhitworthI 
(|2005|), who found complex behaviour for spheroids in the 
non-linear regime. However, the aim of this work is not to 
give an accurate description of the fragment evolution, but 
to obtain an analytical formula which quickly provides in- 
formation about the instability of a given wavelength on the 
shell and which takes into account the effects of external 
pressure. 
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Figure 7. The perturbation growth rate in hydrodynamic simulations, compared to both the thin-shell and the PAGI disper- 
sion relations. The top panels show the simulation with P^xt ~ 10^^^ dynecm"^; the middle panels show the simulation with 



10- 



' dyne cm ^ ; and the bottom panels show the simulation with P^, 



5x10 ^^ dyne cm ^ . In all panels only the 



first (approximately parabolic) peak at low wavenumbers should be compared with the analysis. The other peaks arc higher harmonics, 
which reflect the shapes adopted by fragments in the non-linear condensation regime. The horizontal dashed line at oj = represents the 
separation between unstable and damping modes. 
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Figure 8. Normalised surface density (S/Sq) for the three hydrodynamic simulations with different external pressures: P^ 



10" 



(left), P^ 



10- 



(middle) and P^ 



: 5 X 10 ^'^ dyne cm ^ (right). Panels show a randomly chosen part of the shell surface of the 



size 0.4 X 0.4 radians (the axes are denoted in Hammer projection coordinates). The times are the same as for the middle column of the 
figure[7l i.e. 15, 10 and 4.5 Myr for the low, medium and high pressure simulation, respectively. 
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6 CONCLUSIONS 

We have studied an isothermal ballistic shell, confined from 
both sides by a hot highly rarefied gas having non-zero 
pressure, -P^xt- We have solved the equation of hydrostatic 
equilibrium in a frame of reference moving with the shell 
and shown that the resulting shell density profiles are in 
very good agreement with the profiles obtained from three- 
dimensional hydrodynamic simulations. We have also shown 
that the shell thickness measured at half of its peak density 
agrees well with the thickness of a simple uniform-density 
shell model. 

This allows us to model a fragment forming in the shell, 
due to gravitational instability, as a uniform-density oblate 
spheroid. We use this model to derive a dispersion relation 
for pressure assisted gravitational instability (PAGI). This 
dispersion relation (perturbation growth rate as a function 
of initial fragment wavenumber) depends on the pressure of 
the external medium: the higher Pj,xt ' ^^^ larger the maxi- 
mum unstable wavenumber (i.e. the smaller the smallest un- 
stable fragment), and the faster the growth rate for an unsta- 
ble fragment. If P^xt ~^' ^^^ highest unstable wavenumber 
is 0.6 times smaller than predicted by the standard thin- 
shell analysis using the same parameters; if P^xt =oo, the 
highest unstable wavenumber is 2.2 times larger. The PAGI 
dispersion relation gives approximately the same range of 
unstable wavelengths as the thin-shell dispersion relation if 
PgxT — ^cRiT' where Pq^i^. is the critical external pres- 
sure for which self-gravity and external pressure contribute 
equally to confinement of the shell. 

Finally, we have demonstrated that the predictions of 
the PAGI dispersion relation agree rather well with the 
results of three-dimensional hydrodynamic simulations. In 
particular, (i) the PAGI dispersion relation predicts a max- 
imum unstable wavenumber very similar to the simulations 
(but systematically ~ 25% higher); (ii) modulo this offset, 

the increase in ip^^j with Pj,xt predicted by the PAGI dis- 
persion relation is exactly mirrored by the simulations. 
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